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Inspired by recent advancement of low-power ferroelectic-gated memories and transistors, we 
propose a design of ferroelectic-gated nanoplasmonic devices based on graphene sheets clamped 
in ferroelectric crystals. We show that the two-dimensional plasmons in graphene strongly couple 
with the phonon-polaritons in ferroelectrics at terahertz frequencies, leading to characteristic modal 
wavelength of the order of 100-200 nm at only 3-4 THz. By patterning the ferroelectrics into 
different domains, one can produce compact on-chip plasmonic waveguides, which exhibit negligible 
crosstalk even at 50 nm separation distance. Harnessing the memory effect of ferroelectrics, low- 
power electro-optical switching can be achieved on these plasmonic waveguides. 



The emergence of graphene research in the recent years 
has triggered a significant interest in two-dimensional 
plasmonics. The transport behaviors of graphene re- 
veal extremely low ohmic loss and nearly perfect electron- 
hole symmetryJ^ii The charge-carrier density (or the 
Fermi level) can be conveniently adjusted via chemi- 
cal doping and electrostatic gating, which give rise to 
tunable plasmonic oscillations in the terahertz frequency 
regime. In combination with its remarkable character 
of single-atom thickness and the ability of subwavelength 
light confinement, graphene has become a promising plat- 
form for the new-generation nanoplasmonic devices. So 
far the most widely studied graphene plasmonic struc- 
tures are fabricated on silicon dioxide (Si02) substrates 
or suspended in air,^^^'^ where the Si02 and air serve 
as the dielectric claddings for the creation of surface 
plasmon-polaritons (SPPs).^^ In contrast with conven- 
tional dielectrics, ferroelectrics such as lithium niobate 
(LiNbOa) and lithium tantalate (LiTaOa) bear giant per- 
mittivity and birefringence at terahertz frequencies ^i^i^U 
This peculiarity is associated with the coexistence of ter- 
ahertz phonon-polariton modes^^^^^ and static macro- 
scopic polarization^^"— It is thus of particular interest to 
investigate whether the unique properties of ferroelectrics 
can be integrated into nanoplasmonics,^^^^ and espe- 
cially, display a strong coupling with graphene plasmon- 
ics in the same terahertz range. 

In this paper, we propose a design of ferroelectic-gated 
graphene plasmonic devices operating at THz frequen- 
cies. Figure [1] illustrates the building blocks, while more 
complex structures can be made based on the same idea. 
This type of architecture has at least two appealing fea- 
tures that are not present in the common dielectric- 
graphene-dielectric architecture, and may bring fresh in- 
sight into the design of (classical or quantum) optical 
circuits. The first feature originates from the extremely 
large permittivity and comparatively high Q-value near 
the optical-phonon resonances of ferroelectrics at tera- 
hertz frequencies.^^ When the two-dimensional plasmons 
in graphene become coupled with the phonon-polaritons 
in ferroelectrics by exchanging photons across their inter- 
faces, they form the so-called surface plasmon-phonon- 
polaritons (SPPPs) and lead to around 100 nm modal 
wavelength, even if the driving frequency remains at a 
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FIG. 1. Schematics of a ferroelectric-graphene- ferroelectric 
structure with parallel- and antiparallel-poling domains of the 
macroscopic polarization in the ferroelectrics. Parallel poling 
does not change the Fermi level of graphene from (possible) 
chemical doping, while antiparallel poling can induce an extra 
(positive or negative) electrostatic gating effectively. 



few THz. These subwavelength modes are fundamen- 
tally supported by the atomic-level oscillations of ferro- 
electric ions and are limited by the dissipation through 
anharmonic phonon processes. The second fea- 
ture lies in the large macroscopic polarization in fer- 
roelectrics of the order of 10 jiC cm~^, corresponding 
to a large surface bound-charge density of the order of 
10~^^ cm~^J^— The parallel- and antiparallel-poling 
configurations can induce drastically different free-charge 
densities on graphene, which effectively switch off and on 
the Fermi energy by about 1 eV. These features may 
be employed for constructing nanoplasmonic elements 
(waveguides, resonators, antennae, etc.) through the 
electrostatic gating of ferroelectrics. ^^^^^ Thanks to the 
memory effect of ferroelectrics, an intended design can 
sustain itself without a need of constant input bias, and 
can be conveniently refreshed in a low-power manner, 
similar to those operations in the ferroelectric random- 
access memories (FeRAM) and ferroelectric field-effect 
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transistors (FeFET). 23^24^134 

Let us first study the eigen-modes on an infinite 
ferroelectric-graphene-ferroelectric structure. Suppose 
the materials being stacked along the 2;-axis, and in- 
finitely extended in the xy-plane, matching the same co- 
ordinate system indicated in Fig. [TJ The graphene sheet 
is situated at 2: = with a zero thickness yet a nonzero 
two-dimensional conductivity agi^^^ The two ferroelec- 
tric crystals occupy the semi-infinite regions z < and 
z > 0, respectively. We assume the macroscopic po- 
larizations of the ferroelectrics to be aligned with ei- 
ther -\-z OT —z direction, so the optical axes of the crys- 
tals always coincide with the z-axis.^^^^^ The static po- 
larity mainly affects the electron density on graphene 
but has no appreciable impact on the optical proper- 
ties in bulk crystals. ^3^24,33,34 rpj^^ eigen-solutions of the 
entire structure can be labeled by frequency u and in- 
plane wave numbers kx and ky across all the regions. 
Within each region, the electric field E and the magnetic 
field H are linear combinations of plane waves associ- 
ated with an out-of-plane wavenumber kz. For surface- 
wave solutions, kz is an imaginary number. Due to 
anisotropy, we may define the ordinary and extraordi- 
nary two evanescent wavenumbers with respect to the 
optical z-axis, k!^ = —k^ = (/c^ + ~ (^^/c^)eo, 
= -^l = i^l + ^)^o/^e - (cJ^/c^)eo, respectively, eo 
and ee are the ordinary and extraordinary permittivities 
of the ferroelectrics ^^i^ After some boundary treatment, 
we shall find the dispersion relation (in the cgs units), 

2^ = ^, (1) 

which is an anisotropic generalization to the dispersion 
relation in dielectric-graphene-dielectric structures.^ 

In graphene plasmonics, the Fermi energy measured 
from the Dirac point usually ranges from about ±0.05 eV 
to ±1.0 eV, equivalent to an electron or hole concentra- 
tion Tie varying from 1.84 x 10^^ cm~^ to 7.35 x 10^^ cm~^ 
according to the relation = (lEpl/hvp)'^ /ir, where Vp 
is the Fermi velocity taken as 1 x 10^ cm For 
the experiments below 10 THz, the simple Drude formula 
can quite accurately describe the graphene conductivity 
(in the cgs units) ^i^ii"— cTg = ie'^\Ep\/ nh^ {uj ± i7g ) , where 
7g is the relaxation rate. For an ultra-pure sample, the 
usual mobility limitation comes from the electron scat- 
terings with the thermally-excited acoustic phonons in 
graphene and the remote optical phonons from Si02.^^^^ 
According to the literature, we estimate the mobility for 
our structure to be of the order of 1 x 10^ cm^ V~^s~^ at 
room temperature (~ 300 K), and 1 x 10^ cm^ s~^ 
at low temperature ('^ 100 K).^^ The relaxation rate can 
be calculated from mobility via 7g = ev^/ ii\E^\}^''^^ 

We choose LiNbOa as an example of ferroelectric ma- 
terials. It is especially convenient for low-THz experi- 
ments, because the extraordinarily polarized THz waves 
can be triggered by 800 nm femtosecond laser pulses 
via nonlinear optical response inside LiNbOa.^^'^^ Ac- 
cording to the Raman scattering data,^^^— ^ LiNbOs 



has two fundamental transverse optical-phonon frequen- 
cies cJto and Co^Te, corresponding to the ordinary and 
extraordinary waves, respectively. The main behaviors 
of the permittivities eo and in the 0.1-10 THz fre- 
quency range can be fitted by the Lorentz model, ^^j^ 
eo{(jj) = 60(00) ± (eo(0) - 60(00)) ujI^/ {ul^ - ij? - i7oCj) , 
ee(cj) = ee(oo) ± (ee(0) - ee(oo)) bj\J (cj^g -u? - i7ea;) , 
where 60(00), eo(0), ee(oo), ee(0) are the high-frequency 
and low- frequency limits of eo and ee, 7o and 7e are 
the relaxation rates associated with anharmonic optical- 
phonon decaying. We adopt their room-temperature (~ 
300 K) values from Ref.^^, and set the low-temperature 
(~ 100 K) 7o and 7e to be about one third of their 
room-temperature values according to Ref.^'^ 39 
Co undergo near divergence and sign change around uj^o 
and co^Te, signifying the high-Q resonant coupling be- 
tween photons and optical phonons. They turn back 
to positive at t he longitudin al optical-phonon frequen- 
cies cJlo = ^To\/eo(0)/eo(oo) = 27r x 6.7 THz and uoi^e = 
c^Te\/ee(0)/ee(oo) = 27r X 12.3 THz in this model J^*^ 

For the plane- wave study, we set ky = and choose kx 
as the progressive wavenumber, which can be explicitly 
solved as 

79 / eo(co')ee(co') 9 / . evl \^ 

The modal wavelength A = 27r/Re[/Ca^], the attenua- 
tion length ^ = l/Im[A:a:], and the confining length 
/ = l/Re[A>:e] can be defined accordingly. For a Fermi en- 
ergy between 0.05-1.0 eV, and driving frequency in the 
range of 1-10 THz, the dispersion relation is primarily 
controlled by the second term of Eq. (|2]). The plasmon 
behavior is strongly affected by the optical-phonon reso- 
nances. So the combined excitations are surface plasmon- 
phonon-polaritons (SPPPS). Figure [2] shows the disper- 
sion relation uj{Re[kx]) under different conditions. For 
each given Ep^ the dispersion curves mainly stay in the 
three allowed-bands: uj < cJto, ^lo < oj < c^Te, and 
uj > c^Le, but weakly leak into the two forbidden-bands: 
uJto < oj < cJlo and c^Te < oj < uJi^q. Finite- valued relax- 
ation rates broaden the sharp peaks around the optical- 
phonon resonances. For a small E^^ the dispersion curves 
bend more considerably towards the optical-phonon lines 
and can give very large kx at low frequencies; for a large 
Ep , the dispersion curves mostly attach to the THz light 
line (very close to the cj-axis and cannot be identified 
in the scale of Fig. [2]). In reality, there exist other 
optical-phonon resonances higher than cJto and c^Te,^^^^' 
which will make the curves in Fig. [2] more kinked than 
as shown. But we will only focus on the frequency re- 
gion close to the fundamental optical-phonon resonance 
cJto = 27r X 4.6 THz from below. In Table [U we list 
the calculated characteristic quantities at about 100 K. 
For low-power sensitive THz-photon manipulation, a low 
temperature is helpful for suppressing the thermal noise 
or dissipation. ^-ii^i^i — In Table [H one can see clearly 
the huge effective refractive index and ultra-short modal 
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FIG. 2. Dispersion relations of surface plasmon-phonon- 
polaritons for different Fermi energies on a graphene sheet 
clamped in LiNbOs crystals in the low-temperature 100 K) 
and room-temperature 300 K) conditions. 



TABLE I. Calculated effective refractive index n, modal wave- 
length A, attenuation length ^, and confining length / of 
the surface plasmon-phonon-polaritons with several repre- 
sentative Fermi energies and frequencies at low temperature 
(- 100 K). 
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wavelength compared with the free-space wavelength of 
low-THz photons. Large dissipation occurs at frequen- 
cies above 4.0 THz. But just below it, for ^0.1 eV, 
the wavelength can indeed be squeezed to 100-200 nm 
while the confining length is only about 10-20 nm. 

As an example, we now employ the large difference in 
the length scale of SPPPs under different Fermi energies 



to make low-power subwavelength waveguides. LiNbOs 
is known to possess very large spontaneous polariza- 
tion Pg- In a bulk crystal under zero electric field, 
Pg ~ 70 /iC cm~^,^^'^^'^^ which is equivalent to a sur- 
face bound-charge density ng ~ 4.4 x 10^^ cm~^. In a 
thin film of about 200 nm thick, Pg ^ 5 /iC cm~^^ 
and the equivalent surface bound-charge density is ng ^ 
3.1 X 10^^ cm~^. For our studies, we take the 200 nm 
thickness for each slab of LiNbOs. As shown in Fig. [U 
for the parallel-poling configuration, the bound charges 
of opposite signs from the lower and upper slabs cancel 
each other, leaving an approximately zero-potential set- 
ting for the graphene sheet. Thus the charge-carrier den- 
sity on graphene is solely determined by chemical doping. 
But for the antiparallel-pohng configuration, the bound 
charges of the same sign from the both slabs cause a 
net positive or negative potential on the graphene sheet. 
For the 200 nm thin film clamping, the induced charge- 
carrier density is Uc ~ 2ng = 6.2 x 10^^ cm~^ which is 
equivalent to a nearly ±1 eV electrostatic gating based 
on the preceding calculation. In the numerical simu- 
lation below, we assume that a small Ep = 0.05 eV 
is built in graphene from chemical doping, and is re- 
served between the parallel-poling domains, while a large 
Ep ^ 0.05 ± 1.0 ^ ±1.0 eV is generated between the 
antiparallel-pohng domains due to the ferroelectric gat- 
ing. The slight difference in \Ep\ between the electron 
and hole cases is neglected because of the smallness of 
0.05 eV compared with 1.0 eV. The electron- hole sym- 
metry in graphene plays an important role here. 

We choose the frequency 3.5 THz for all our finite- 
difference time-domain (FDTD) simulations. The free- 
space photon wavelength is 85.7 /im. For the domain pol- 
ing configuration shown in Fig. [3l a single E^ = 0.05 eV 
channel is produced in between two Ep = 1.0 eV barri- 
ers. As can be inferred from our previous discussion, a 
low-£^F channel is more "dielectric" in the sense that it 
carries less charges but hosts more photons, whereas a 
high-£^F channel is more "metallic" in the sense that it 
carries more charges but permits less photons. Thus the 
deep subwavelength SPPP modes preferably flow in the 
middle channel with very tiny penetration into the left- 
and right-barriers. The lower panel of Fig. [3] shows the 
real and imaginary parts of the effective refractive index 
h changing with the channel width. One can see that the 
waveguide mode undergoes exponentially stronger atten- 
uation and weaker subwavelength after the channel width 
goes down to below 100 nm, consistent with the 93 nm 
plane- wave modal wavelength at 3.5 THz in Table [H We 
choose 100 nm to plot the profiles of each electric-field 
component in the upper panel of Fig. [3l For the effec- 
tive refractive index n = 810.81 + 43.10i in this case, the 
modal wavelength A = 105.6 nm, and the attenuation 
length ^ = 316.2 nm. 

The exceptional confining quality of the prescribed 
ferroelectric-gated channel can be manifested by putting 
two such channels close to each other and see how the 
modes in the two channels become coupled by varying 
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n= 810.81 + 43. lOi 
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FIG. 3. Simulated effective refractive index and electric- field 
profile of the fundamental waveguide mode in a parallel-poling 
channel between two antiparallel-poling barriers. (Blue and 
red arrows indicate domain orientations.) 



the separation distance ^iii^i^^ We may define a dimen- 
sionless crosstalk number, 



C = -101o£ 



10 



■ Re \h. 



4(Re[na] -Re[ns]) 



(3) 



which estimates (in terms of dB) the number of modal 
(not free-space) wavelength needed, for an injected power 
initially in one channel to be transferred into the other 
and then be transferred back. In Fig. 21 we can see the 
simulated crosstalk number to be mostly in the range 
of —45 to —50 dB. For the plotted pattern of 20 nm 
separation distance, this number is still far below zero 
at —29.5 dB. The dissipation has certainly come in 
at a much shorter propagation distance, so the waves 
cannot really travel that far.^^ But this number still 
shows the extreme confining quality in these ferroelectric- 
graphene waveguides compared with conventional dielec- 
tric and plasmonic waveguides. For example, silicon 
waveguides with similar dimensions at infrared frequen- 
cies have a crosstalk number of the order of —10 dB in 
accordance with our definition. Bao.sSro.sTiOa-metal 
inter-layer plasmonic waveguides at visible-light frequen- 



cies can have a crosstalk number of the order of —30 dB, 
but must take a much larger separation distance. One 
may notice a singular drop at about 50 nm on the cal- 
culated curve in the lower panel of Fig. SI Based on an 

h 1200 nm -| 
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FIG. 4. Simulated crosstalk number and electric-field profiles 
of the antisymmetric and symmetric waveguide modes in two 
100 nm wide parallel-poling channels. (Blue and red arrows 
indicate domain orientations.) 



analysis to the dipole-dipole coupling in this particular 
system, we find it to be due to a competition between 
the relative magnitudes of field components. A large Ey 
field tends to make the symmetric mode have a higher 
refractive index, while a large Ex or Ez field tends to 
make the antisymmetric mode have a higher refractive in- 
dex. The 50 nm separation distance between two 100 nm 
wide waveguides happens to be the turning point, where 
Re[na] — Ke[ns] — 0. A more thorough study on this 
phenomena will be performed. 

We acknowledge the financial support by NSF (EGGS 
Award No. 1028568) and the AFOSR MURI (Award No. 
FA9550-12-1-0488). 



* |nicfang@mit.edu| 



^ M. Jablan, H. Buljan, and M. Soljacic, Phys. Rev. B 80, 



5 



245435 (2009). 

^ J. Long, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, 
H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, and 
F. Wang, Nature Nanotechnology 6, 630 (2011). 

^ S. Ryu, J. Maultzsch, M. Y. Han, P. Kim, and L. E. Brus, 
ACS Nano 5, 4123 (2011). 

^ F. H. L. Koppens, D. E. Chang, and F. J. Garcia de Abajo, 
Nano Lett. 11, 3370 (2011). 

^ A. R. Davoyan, V. V. Popov, and S. A. Nikitov, Phys. 
Rev. Lett. 108, 127401 (2012). 

^ Q. Bao and K. P. Loh, ACS Nano 6, 3677 (2012). 

^ A. Vakil and N. Engheta, Science 332, 1291 (2011). 

^ J. Chen, M. Badioh, P. Alonso- Gonzalez, S. Thongrat- 
tanasiri, F. Huth, J. Osmond, M. Spasenovic, A. Cen- 
teno, A. Pesquera, P. Godignon, A. Z. Elorza, N. Camara, 

F. Javier Garcia de Abajo, R. Hillenbrand, and F. H. L. 
Koppens, Nature 487, 77 (2012). 

^ Z. Fei, A. S. Rodin, G. O. Andreev, W. Bao, A. S. 
McLeod, M. Wagner, L. M. Zhang, Z. Zhao, M. Thiemens, 

G. Dominguez, M. M. Fogler, A. H. C. Neto, C. N. Lau, 
F. Keilmann, and D. N. Basov, Nature 487, 82 (2012). 
A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

K. L Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, and 
P. Kim, Phys. Rev. Lett. 101, 096802 (2008). 
S. V. Morozov, K. S. Novoselov, M. L Katsnelson, 
F. Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, 
Phys. Rev. Lett. 100, 016602 (2008). 
^2 N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 

5. Das Sarma, S. Adam, E. H. Hwang, and E. Rossi, Rev. 
Mod. Phys. 83, 407 (2011). 

J. H. Chen, C. Jang, S. Xiao, M. Ishigami, and M. S. 

Fuhrer, Nature Nanotechnology 3, 206 (2008). 

A. Y. Nikitin, F. Guinea, F. J. Garcia- Vidal, and 

L. Martm-Moreno, Phys. Rev. B 84, 161407 (2011). 

J. Christensen, A. Manjavacas, S. Thongrattanasiri, 

F. H. L. Koppens, and F. J. Garcia de Abajo, ACS Nano 

6, 431 (2012). 

E. N. Economou, Phys. Rev. 182, 539 (1969). 
T. Feurer, N. S. Stoyanov, D. W. Ward, J. C. Vaughan, 
E. R. Statz, and K. A. Nelson, Annu. Rev. Mater. Res. 
37, 317 (2007). 

^° Y. M. Sun, Z. L. Mao, B. H. Hou, G. Q. Liu, and L. Wang, 
Chin. Phys. Lett. 24, 414 (2007). 

N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Brooks Cole, 1976). 

V. Marek, First-principles study of ferroelectric oxides^ 
Ph.D. thesis, Universite de Liege (2003). 



N. Setter, D. Damjanovic, L. Eng, G. Fox, S. Gevorgian, 
S. Hong, A. Kingon, H. Kohlstedt, N. Y. Park, G. B. 
Stephenson, I. Stolitchnov, A. K. Taganstev, D. V. Taylor, 
T. Yamada, and S. Streiffer, J. Appl. Phys. 100, 051606 
(2006). 

M. Dawber, K. M. Rabe, and J. F. Scott, Rev. Mod. Phys. 
77, 1083 (2005). 
2^ S. Sanna and W. G. Schmidt, Phys. Rev. B 81, 214116 
(2010). 

2^ S. W. Liu and M. Xiao, Appl. Phys. Lett. 88, 143512 
(2006). 

J. E. Spanier, A. M. Kolpak, J. J. Urban, L Grinberg, 
L. Ouyang, W. S. Yun, A. M. Rappe, and H. Park, Nano 
Lett. 6, 735 (2006). 

M. J. Dicken, L. A. Sweatlock, D. Pacifici, H. J. Lezec, 
K. Bhattacharya, and H. A. Atwater, Nano Lett. 8, 4048 
(2008). 

A. Harhira, L. Guilbert, P. Bourson, and H. Rinnert, Phys. 
Status Solidi C 4, 926 (2007). 
^° V. H. Arakehan and N. M. Hovsepian, Phys. Status Solidi 
B 164, 147 (1991). 

A. Haussmann, P. Milde, C. Erler, and L. M. Eng, Nano 
Lett. 9, 763 (2009). 

D. Li and D. A. Bonnell, Annu. Rev. Mater. Res. 38, 351 

(2008) . 

Y. Zheng, G.-X. Ni, C.-T. Toh, M.-G. Zeng, S.-T. Chen, 
K. Yao, and B. Ozyilmaz, Appl. Phys. Lett. 94, 163505 

(2009) . 

E. B. Song, B. Lian, S. M. Kim, S. Lee, T.-K. Chung, 
M. Wang, C. Zeng, G. Xu, K. Wong, Y. Zhou, H. L Rasool, 
D. H. Seo, C. H.-J, J. Heo, S. Seo, and K. L. Wang, Appl. 
Phys. Lett. 99, 042109 (2011). 

X. Hong, A. Posadas, K. Zou, C. H. Ahn, and J. Zhu, 
Phys. Rev. Lett. 102, 136808 (2009). 
J. Hebling, K.-L. Yeh, M. C. Hoffmann, and K. A. Nelson, 
IEEE Journal of Selected Topics in Quantum Electronics 
14, 345 (2008). 

S. V. Ivanova, V. S. Gorelik, and B. A. Strukov, Ferro- 
electrics 21, 563 (1978). 

P. Capek, G. Stone, V. Dierolf, C. Althouse, and 
V. Gopolan, Phys. Status Solidi C 4, 830 (2007). 
M. L. Hua, C. T. Chia, J. Y. Changa, W. S. Tse, and 
J. T. Yu, Mater. Chem. Phys. 78, 358 (2002). 
^° V. Joshi, D. Roy, and M. L. Mecartney, Appl. Phys. Lett. 
63, 1331 (1993). 

C. Yeh and F. Shimabukuro, The Essence of Dielectric 
Waveguides (Springer, 2008). 

G. Veronis and S. Fan, Opt. Express 16, 2129 (2008). 



